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This paper addresses the problem of vortex formation during a rapid quench in a superconducting 
film. It builds on previous work showing that in a local gauge theory there are two distinct mecha- 
nisms of defect formation, based on fluctuations of the scalar and gauge fields, respectively. We show 
how vortex formation in a thin film differs from the fully two-dimensional case, on which most the- 
oretical studies have focused. We discuss ways of testing theoretical predictions in superconductor 
experiments and analyse the results of recent experiments in this light. 



I. INTRODUCTION 

There has been much recent interest in testing cosmo- 
logical defect-formation scenarios 0, 0> Q in the labo- 
ratory. In the Nineties, experiments in liquid crystals 
ULIMlI an d especially in liquid 3 He @,H1 clearly showed 
the formation of defects in line with predictions. In 4 He, 
experiments have proved negative but there is good 
reason to believe that this is due to the fact that vortices 
disappear too rapidly to be seen 0, ^| . 

More recently, the focus of the ex peri mental acti vity 
has shifted towards superconductors |l2t fl3l PiU ITU ITa . 
11^. Il8j . In the pioneering experiment by Carmi and 
Polturak |l2j | , no evidence of spontaneous vortex forma- 
tion was seen, but later experiments with a multiple- 
Josephson-j unction l oop [lflj and with an annular Joseph- 
son junction 0, 0, 0] have been positive. Earlier this 
year, successful experiments have also been carried out 
with an array of superconducting rings |17| and a thin 
superconducting film [Isj . 

It is not possible to apply the original estimates of de- 
fect numbers directly to the case of the superconducting 
transition, because it involves a local gau ge theory rather 
than a global symmetry breaking 0, [H| |2(J . In ref . ^3 > 
a somewhat ad hoc modification of the standard predic- 
tion was proposed for comparison with experiment, but 
a more firmly based estimate is needed. The problem of 
defect formation in a gauge theory has been analyzed by 
Hindmarsh and Rajantie 3], who showed that in addition 
to the standard mechanism of defect formation, there is 
a second mechanism operating. This flux trapping mech- 
anism stems from the fluctuations in the magnetic field, 
which get frozen out during the quench. 

The two mechanisms produce defect distributions with 
characteristically very different correlation properties. In 
a gauge theory, one would expect positive vortex-vortex 
correlations at short distances. In other words, the vor- 
tices should be formed in clusters of equal sign. This 
is supported by numerical simulations carried out in the 
Abelian Higgs model pit 12^ . On the other hand, the 



vortices should have negative correlations at all distances 
in global theories. 

One of us has studied an SO (3) gauge theory, 
and shown that monopolcs form, even when there is no 
true phase transition, only a smooth crossover. The 
monopole distribution again shows positive monopole- 
monopole correlation, in contrast to the situation in 
a global symmetry theory, where one would expect 
monopole-antimonopole correlations. 

In this paper, we examine the spontaneous formation 
of vortices (or fluxons or fiuxoids) in a thin film undergo- 
ing the normal-to-superconducting transition in the light 
of existing and future experiments. This setup differs 
significantly from the fully two-dimensional case studied 
in Ref. [2^. because the magnetic field extends outside 
the film and leads to long-range 1 jr interaction between 
vortices |2 lj. Therefore, there are actually many simi- 
larities with monopole formation in three dimensions. In 
both cases, the interaction is screened by thermal effects, 
although the screening is less effective here. 

We also discuss the results of the experiments carried 
out so far, and conclude that whilst they are compatible 
with the two mechanisms, they do not currently consti- 
tute a positive confirmation of either. 

Throughout this paper, we use the natural units c = 
% = he = 1- We will write three- vectors with bold-faced 
symbols as k, and two-vectors on the film with an arrow 
as k. 



II. TWO-DIMENSIONAL 
SUPERCONDUCTORS 

Theoretical studies of defect formation in two dimen- 
sions have shown that there are essentially two different 
regimes. At a weak enough coupling or low enough crit- 
ical temperature, the relevant degree of freedom is the 
order parameter field and the process is well described by 
the Kibble-Zurek mechanism [HQ. On the other hand, 
at strong coupling or high critical temperature, it is dom- 
inated by thermal fluctuations of the gauge field . 
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The simplest example of the Kibble-Zurek mechanism 
is the phase transition in a Ginzburg-Landau model with 
a global U(l) symmetry. Because the mechanism is very 
generic, we can consider a two-dimensional relativistic 
version of the model, with the Lagrangian density 
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The U(l) symmetry, which is spontaneously broken at 
zero temperature, is restored at temperatures above a 
critical value T c , Above T C) the field has a finite correla- 
tion length, which diverges at T c . In the broken phase, 
the phase angle is a Goldstone mode and therefore it has 
an infinite correlation length. 

Imagine now that the temperature is initially above T c , 
and is gradually decreased so that the system undergoes 
a phase transition. When T c is approached, the correla- 
tion length would have to grow faster and faster in or- 
der to stay in equilibrium. Sooner or later that becomes 
impossible. As a result, the phase angle can only be or- 
dered over some finite distance £ in the broken phase, 
and this leads to formation of vortices with number den- 
sity n ~ l/£ 2 pj. The length scale £ has a power-law 
dependence on the cooling rate, with an exponent that 
depends on the universality class of the transition [2| . 

The spatial vortex distribution can be characterized 
by how the typical net vortex number AA inside a cir- 
cular curve depends on its radius R. As long as R 3> £, 
the curve should pass through roughly 2irR/£ correlated 
domains, and each time it moves from one domain to 
another, the phase changes by a random amount. As a 
consequence, the net vortex number should behave like a 
random walk of 2nR/^ steps, 
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For comparison, if the signs of the vortices were ran- 
domly distributed, the net vortex number would grow 
like AN(R) oc R. This shows that the vortices formed 
by this mechanism have negative correlations. 

Superconductors have a spontaneously broken local 
gauge symmetry, and therefore we modify the Lagrangian 
by adding a gauge field, i.e., the electromagnetic vector 
potential A, 
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C = -[E Z - B 2 ) +<tf<l>-D<f-D<t>-l3(<l>'<l>-v'y, (3) 

where we have introduced the covariant derivative D = 
V + ieA, and used the gauge invariance to set the electric 
scalar potential to zero. With this choice, the electric 
and magnetic fields are given by E = —dA/dt and B = 
d 1 A 2 -d 2 A 1 . 

This describes a relativistic, fully two-dimensional su- 
perconductor, and is not a good approximation for a su- 
perconductor film, as we will see later. Nevertheless, we 
will use it to illustrate the flux trapping mechanism of 
defect formation Q. 



Classically, the initial thermal state is given by the 
Boltzmann distribution p[(j), AJ oc exp(— E[<j), A^j/T), 
where 
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The most relevant terms for us are the one containing 
the magnetic field and the covariant derivative term. 

If the magnetic field vanishes, the covariant derivative 
reduces to the ordinary gradient, and one expects that 
the above Kibble-Zurek picture is valid. However, at a 
non-zero temperature, there are always thermal fluctua- 
tions, which have to be taken into account. In the sym- 
metric phase, the effect of the scalar <j> on the magnetic 
field is insignificant, and to a good approximation we can 
say that B is a Gaussian random field with a two-point 
function 



or 



(B(x)B(y))nT6(x-y), 
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This is simply the classical Ray leigh- Jeans spectrum of 
blackbody radiation, and therefore it is only valid when 
\k\ < k B T/ch S3 440(r/K)m" 1 . This constraint may be 
important in the case of low-temperature superconduc- 
tors, but we will not discuss it here any further.. 

In the broken phase, the Meissner effect suppresses the 
fluctuations, and the two-point function (for modes with 
k<T) becomes 



(B(k)B(q)) S3 T 



k 2 + A 



— (27T) 2 S(k + q). 
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where A is the penetration depth, which starts to decrease 
from infinity at T c . 

However, the magnetic field is conserved, because 
d Q B = V • E, where E = {-E 2 ,E 1 ). Therefore, the 
long-wavelength modes can only decay very slowly. Be- 
cause of causality, the decay rate must be lower than the 
wave number "f(k) < k, but there are also other, much 
stronger constraints arising from the detailed dynamics. 

In any case, for any cooling rate, there are going to be 
some long-wavelength modes that are too slow to react. 
Their amplitudes are therefore frozen to the values they 
had before the transition. The corresponding magnetic 
flux gets trapped into vortices. We can define a critical 
wavenumber k c , so that this happens to modes with k < 
k c . 

Now, consider the typical (root-mean-squared) mag- 
netic flux through a circle of radius R S3 l/k c , 
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The whole contribution to this quantity comes from 
modes with k < k c , and therefore Q(R) does not change 
during the transition. After the transition, it will have 
the same value as in the symmetric phase before the tran- 
sition, and we can calculate that, using Eq. @ (and ne- 
glecting factors of order 1), 



$(R) « Vr/k c . 



(9) 



However, because of Meissner effect, this magnetic flux 
must be confined into vortices. Each vortex carries one 
flux quantum <E>o = 2ir/e, and therefore we can say that 
the typical number of vortices in a region of radius R = 
l/k c is 



N(R) = e$(R) fa eVT/k c . 



(10) 



Dividing this by the area of such a region gives an esti- 
mate for the number density Q 



zkrVT. 



(11) 



Nevertheless, the most obvious difference with the Kib- 
ble mechanism is that all N(R) vortices formed in any 
given region have the same sign, which leads to positive 
correlations between vortices and formation of clusters of 
N(R) vortices. These clusters have been seen in numer- 
ical simulations in Ref. [22| . 

Note that our discussion was only based on the be- 
haviour of the long- wavelength magnetic fields, namely 
the qualitative form of the two-point function in Eq. (JJ, 
causality and the conservation of magnetic field, ft 
should therefore be independent of the microscopic de- 
tails of the system and should apply to real superconduc- 
tors as well as to the toy model in Eq. J2J. However, as 
we shall see, one will have to consider a three-dimensional 
setup. 



III. SUPERCONDUCTING FILMS 

Superconductor experiments are typically carried out 
with a thin superconducting film, but the fully two- 
dimensional geometry discussed in Section [H] does not 
describe that situation correctly, because the magnetic 
field should extend outside the superconductor and in- 
duce long-range interactions between vortices |24| . 

In the symmetric phase, we can calculate the analogue 
of Eq. Q, i.e., the two-point function on the film from the 
three-dimensional correlation function in the Feynman 
gauge, 



<A J (k)^(q)) = ^(2 7 r) 3 ( 5(k + q 
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Again, this is equivalent to the classical Rayleigh-Jeans 
spectrum. For the magnetic field normal to the film, we 
find 



(B 3 (x,0)B 3 (y,0))=T 



d 3 k P 
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,ik-{x-y) 



(13) 



where we have set k = (k, ka). Equivalently, defining 

B{k) = [ d 2 xB 3 (x,0)e-^ s , (14) 



we have 



Tk 

(B(k)B(q)) = —(2ir) 2 5(k + q). 



(15) 



It is less straightforward to calculate the two-point func- 
tion in the broken phase, because the gauge field is mas- 
sive on the film but massless outside it. Therefore, we 
use a different approach. Deep in the broken phase, the 
vortices are well defined, localized objects, and we can ap- 
proximate them by point particles on the film. Because 
the magnetic field cannot penetrate the film, we can flip 
the direction of the magnetic field on one side of the 
film, without changing the energy of the configuration. 
This changes a vortex into a magnetic monopole of charge 
<7m = 2<E»o = 4ir/e. Therefore, the vortices must behave 
as magnetic monopoles confined on the film. Correspond- 
ingly, their interaction potential falls off as 1/r, as op- 
posed to the log(r) of global vortices and exp(mr) of local 
vortices in a fully two-dimensional case. In fact, this re- 
sult is only valid at distances longer than the Pearl length 
Ap« \ 2 /d, where d is the thickness of the film [24|. In 
what follows, we assume that Ap is shorter than any of 
the length scales we discuss. 

We now assume that the number density of thermal 
vortices is nt = n + + ri—, where n + and n_ are the den- 
sities of positive and negative vortices, respectively, and 
calculate what kind of correlations the magnetic Coulomb 
interaction induces. This calculation is very similar to 
the calculation of Debye screening in three dimensions. 
In the absence of thermal fluctuations, the magnetic po- 
tential around a vortex would be simply cf>Q (r) — qu /Airr. 
However, the presence of such a potential biases the ther- 
mal distribution of vortices and antivortices such that 



n±(r) 



fexp( T ^ W /T)^ (l T ^). 

(16) 



The magnetic flux density is 
B(r) = q M (n + (r) - ra_(r)) = - 



T 



A, 



where we have defined 



A 



T 
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This flux density, in turn, modifies the potential, and we 
find the integral equation 

0(f) = qM | f fr^E?) - J_ [ffM?l 
47rr J 4:7r\r-r'\ 47rr A scr J ^\r-r'\'' 

' : ' !) ' 

We can solve this equation in Fourier space, and find 



1 gM 
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(20) 
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In coordinate space, 4> behaves asymptotically as 
<t>{r) ~ t^- at r < A scr , 



(r) 



2 

scr 



4 7rr 3 



at 



> A s 



which shows that the thermal fluctuations screen the in- 
teraction potential, but the potential still has a power-law 
form rather than an exponential. Including the contribu- 
tion from the original vortex, the magnetic flux density 
around a vortex is 



A, 



(21) 



and since there is a density rif of these vortices each with 
charge qM, the two-point function of the flux density is 



(B(x)B(y)) =qn f B(\x-y\). 
In Fourier space, this corresponds to 



(22) 



(B(k)B(q))=qn f 
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2 fc + A s 



(27T) 2 5(k + q) 
-{2Tr) 2 S(k + q). (23) 



We can rewrite this as 



(B(k)B(q)) 



T 



2 1 + fcA, 



(2ir) 2 5(k + q), 



(24) 



where we have used Eq. JTSJ. From this expression, it 
is easy to see that the long-distance (low k) behaviour 
agrees with the symmetric phase, where A scr is small. 

The crucial difference between Eqs. {7} and l|24|) is that 
in the case of a film, the important length scale is not the 
penetration depth A (or the Pearl length Ap), but the 
screening length A scr . Whilst A starts from infinity and 
approaches the zero-temperature penetration depth, A scr 
behaves in the opposite way: It is initially very short and 
approaches infinity at zero temperature. 

As A scr increases, each mode gets suppressed only 
when the wavelength becomes shorter than the screen- 
ing length, l/k<A scr . Therefore, the long-wavelength 
modes feel the effect of the phase transition later than 
short-wavelength modes. In fact, the present case is 
actually very similar to monopolc formation in three- 
dimensional non-Abelian theories |23| , since the magnetic 
charge-charge correlator in that case, 



(/J M (k)p A /(q)> wT 



l + k 2 /M- 



■(27r) 3 <5(k- 



(25) 



behaves in the same way. 

In any case, Eq. I|24fl shows that if the system could 
stay in equilibrium, the fluctuations of all wavelengths 
would be suppressed in the zero-temperature limit, where 



A scr — > co. However, the slowness of the long- wavelength 
modes implies that the modes with wave numbers below a 
certain critical value k c freeze out to the amplitudes they 
had in the symmetric phase, and we can approximate 

Tk 

(B(k)B(q)) = — exp(-fc/fc c )(27r) 2 (5(fc + q). (26) 

At this point, we do not specify the value of k c , but only 
note that it must be finite and non-zero. The precise form 
of the cutoff at k c is irrelevant, and this exponential form 
simplifies the calculations, because the coordinate space 
correlator has the simple form 



(B(0)B(r)) 



Tkl 



4tt (l + r 2 fc2)5/2- 



This is constant at short distances, r -C k c 
that we can deduce the number density 



(27) 



and from 



n^^BPW)^^. (28) 

This can also be seen by calculating the variance of the 
flux <&(-/?) through a circle of radius R, 



(* 2 > 



d 2 x / d 2 y(B(x)B(y)) 



{TT/2)Tk 3 c R 4 , 
TRln k c R, 



for R < k t 
for R > k, 



(29) 



At short distances, the typical (rms) flux through a small 
curve is proportional to its area and implies a uniform 
flux density B w (Tk 3 ) 1 / 2 , in agreement with Eq. 
At long distances, the scaling is the same as in the Kibble- 
Zurek scenario, AN(R) ~ i? 1 / 2 , apart of the logarithm. 

By choosing R = k^ 1 in Eq. (|2l?|> we can also see that 
the number of vortices in a cluster is 



(30) 



These results agree with the estimates in Ref. || . 

So far, we have not discussed the value of k c , and in- 
deed, it is generally more sensitive to the detailed dynam- 
ics than the above considerations. On the other hand, 
one can carry out an A: c -independent test of the mecha- 
nism by measuring, for instance, the combination N^n, 




which is predicted to be 



e 4 T 2 
(4tt) 3 



(31) 



To get a crude estimate of k c , let us assume that Ohm's 
law j — oE is valid, and that the dynamics of the mag- 
netic field fluctuations is dominated by the conductivity 
(j. In that case, Maxwell's equations imply that the dy- 
namics of the magnetic field normal to the film is given 
by the damped wave equation 



d 2 B z 
dt 2 



d 2 B z + V 2 B Z 



dBz 
dt 



(32) 
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Ignoring the z derivative, which would give rise to ra- 
diation, we find that the maximum decay rate of any 
long- wavelength mode with wavenumber |fc| <C a is 
7(fc) = k 2 /a. According to Eq. I|24|l . the amplitude of 
the mode must decrease as 

B(k) ~ ^T/A scr (33) 

for the system to stay in equilibrium once A scr > 1 jk. The 

time scale for decrease of B{k) is therefore the quench 
time tq. The decay rate 1/tq exceeds the maximum 
decay rate 7(fc), if 

k < K « yfij7 Q . (34) 

IV. COMPARISON WITH EXPERIMENT 

We now turn to the question of whether these results 
can throw any light on recent experimental results. Ex- 
periments that were very similar to the setup we have dis- 
cussed have been carried out by Carmi and Polturak ^3] 
and more recently by Maniv et al [Tsj]. They involved a 
rapid quench of a thin film of the high-temperature su- 
perconductor YBa 2 Cu307_5, using a SQUID to measure 
the net flux following the quench. This corresponds to 
the net number of vortices minus antivortices produced 
by the quench. 

In the first set of experiments |12| . no evidence for 
vortices were found. The standard causality argument 
used in cases of global symmetry breaking 0, sug- 
gests that the expected scalar-field correlation length £ 
at the relevant time after the transition should in this 
case be about 0.1 /im. In the 1cm 2 area of the sample, 
the number of vortices N should then be about 10 9 or 
10 10 . If one assumes, with Zurek 2], that the phase of 
the scalar field makes a random walk around the perime- 
ter of the sample with a step length of £, then one finds 
AN ~ A 1 / 4 , suggesting that AN should be at least of 
order 100. The experiment should have seen anything in 
excess of A A ~ 20, so this was thought to be a contra- 
diction. 

In the revised experiment ^8|j faster cooling rates of 
around dT/dt ss — 10 8 K/s were used. This increases 
the predicted number of vortices, and indeed, net vortex 
numbers between 30 and 45 were measured. This num- 
ber showed a weak dependence on the cooling rate, with 
faster quenches producing more vortices. This is com- 
patible with theoretical predictions. The Kibble-Zurek 
scenario with overdamped dynamics would yield a do- 
main size of £ oc | dT/dt \ ~ 1 / 4 . The net flux would then 
be 

$ = AJV* « (i?/0 1/2 $ oc \dT/dt\ 1/8 , (35) 

where R sa 1 cm is the linear size of the film. The ex- 
perimental results are consistent with this scaling [l8| . 
although the error bars are so large that one cannot use 



that data to determine the exponent. The prefactor in 
the theoretical prediction is rather uncertain, and to ob- 
tain a fit the authors reduced the prediction by a factor 
of four. Extrapolation to the cooling rate 20 K/s used 
in Ref. |l2| then gives AN 6, which is below the noise 
level and explains why no vortices were seen in the orig- 
inal experiment. 

However, Carmi and Polturak |12| rightly pointed out 
that these arguments are unreliable in the case of a gauge 
theory, and suggested an ad hoc modification leading to a 
prediction AN ~ 0.1N 1 / 2 ~ 10 4 , which is m very severe 
disagreement with the experimental results. 

We wish to argue, however, that this modification is 
not the appropriate way of treating a gauge field. It is 
important to note that the two mechanisms of defect for- 
mation, based on fluctuations of the scalar field and the 
gauge field, are characterized by different length scales. 
In the case of the scalar field, it is the correlation length 
£; for the gauge field, it is l/k c . So to estimate the contri- 
bution to be expected from magnetic-field fluctuations, 
we need to estimate k c . This is hard to do, because it is 
not clear how fast gauge-field correlations can propagate. 
However, we can still put a limit on the expected number 
of defects formed. It seems probable that l/k c would be 
smaller than the size of the sample. In that case, the net 
flux through the sample should be essentially the same 
as that in the high-temperature phase due to thermal 
fluctuations. From Eq. (|29|l . we find 

A A" « Ve 2 TR. (36) 

With a sample size R w 1 cm, this gives AN ~ 10. If 
in fact l/k c is larger than R, then the net flux would 
actually be less than this, 

AN a v/e 2 Tfc c 3 i? 4 oc \dT/dt\ 3/ \ (37) 

where we have used Eq. I)34[l. 

The above estimate shows that the contribution from 
the gauge-field fluctuations to the net flux is so small that 
it would have been below the noise level in Ref. ^2]> J us t 
like the contribution from the Kibblc-Zurck mechanism. 
However, in contrast to the Kibble-Zurek mechanism, the 
contribution is independent of the cooling rate, and the 
faster cooling used in Ref. would not help detecting 
it. At most, it would shift the data points in Fig. 4 of 
Ref. up by a small amount, but because the errors 
are relatively large and the overall scaling of the theoret- 
ical prediction the data points are compared with is un- 
determined, one cannot draw any conclusions about the 
gauge-field contribution based on these measurements. 

We conclude therefore that it is unlikely that 
the magnetic-field fluctuations contribute a detectable 
amount to AN. This does not mean, however, that 
they can be ignored. The two mechanisms differ radi- 
cally in the nature of the resulting defect distributions. If 
scalar- field fluctuations are dominant, one expects strong 
vortex-anti-vortex correlations — the two-dimensional 
analogue of a preponderance of small vortex loops - 
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whereas, as noted above, if magnetic-field fluctuations 
dominate, we expect clusters of similarly oriented defects. 

Clearly, it would be very desirable to be able to get 
some information on the distribution of vortices, and in 
particular to measure N as well as AN. Only in that way 
could one separate the contributions of the two mecha- 
nisms. Achieving this in an experiment with a super- 
conductor film is probably difficult, because the vortices 
produced during the quench may escape the sample in a 
short time. The flux distribution would therefore have to 
be measured very rapidly immediately after the quench. 

However, a recent experiment by Kirtley et al [T?) 
shows how this problem can be solved. Instead of a film, 
they used a square array of 12x12 thin (50 nm) super- 
conducting rings with inner radius of 10 /xm and outer 
radius of 15 /im, spaced by 60 /zm. After the quench, 
the magnetic flux through any ring can only change if a 
vortex can cross the high energy barrier created by the 
superconducting ring, and therefore the original fluxes 
survive a long enough time to be measured with a scan- 
ning SQUID. 

For this experiment, the Kibble-Zurek mechanism 
would predict that each ring behaves independently of 
the others. From Eq. (J2J), the typical number of fluxons 
trapped in a ring would be 

N (x\dT/dt\ 1/8 . (38) 

In fact, there was typically only one flux quantum or no 
flux at all through one ring in the experiment, and it is 
therefore not at all clear one can use Eq. H38fl to predict 
the fluxes. Nevertheless, if one interprets it as the proba- 
bility of having a fluxon, the result of the experiment was 
in a clear disagreement with it. Instead of a power-law, 
the measurements showed a very weak, perhaps logarith- 
mic, dependence on the cooling rate. 

Kirtley et al j^J were able to fit their data with a 



simple model, which is similar in spirit to the flux trap- 
ping scenario. They argue that any vortices formed by 
the Kibble-Zurek mechanism would be 'washed out' by 
their thermally activated vortex mechanism. Because the 
rings are relatively far apart from each other, they can 
be assumed to be independent, and the problem simpli- 
fies considerably. For this reason, the results in Ref. [13 
cannot be thought to have confirmed the flux trapping 
scenario. 

On the other hand, it seems that the predictions for the 
spatial vortex distribution, which can best distinguish the 
two mechanisms, could be tested with a relatively simple 
experimental setup. Instead of a homogeneous film or 
an array of rings, one would use a film with small holes 
in a regular array. The vortices formed during a quench 
would simply move to the nearest hole and get trapped 
there. They could then be measured with a scanning 
SQUID device in the same way as in Ref. |17J . 

It should be noted, however, that there are many other 
complicating factors in the case of a superconducting 
film, notably the importance of impurities and pinning 
centres, which could drastically affect the vortex num- 
bers. Moreover, high-temperature superconductors in- 
volve Cooper pairs with d-wave pairing, so a full descrip- 
tion should involve a multi-component order parameter. 
As a result, the vortices can have more structure, for 
example antiferromagnetic cores. The effects of these 
complications on defect formation remain to be studied. 
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